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Abstract We present a theory for the full distribution of current fluc- 
tuations in incoherent diffusive superconducting junctions, subjected to a 
voltage bias. This theory of full counting statistics of incoherent multiple 
Andreev reflections is valid for arbitrary applied voltage. We present a de- 
tailed discussion of the properties of the first four cumulants as well as the 
low and high voltage regimes of the full counting statistics. The work is an 
extension of the results of Pilgram and the author, Phys. Rev. Lett. 94, 
086806 (2005). 
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1 Introduction 

The mechanism for charge transport in voltage biased superconducting junc- 
tions is Multiple Andreev Reflections (MAR) [T][2]. Coherent MAR-theory 
[3l|4l[5] has provided a quantum mechanical picture for current transport in 
mesoscopic superconducting junctions. The predictions of the theory were 
found to be in impressive agreement with experimental results [SjCfM] • This 
development, with the emphasis on the properties of atomic point contact 
junctions, was recently reviewed [3]. Additional insight in coherent MAR 
transport, as e.g. the effective charge transferred, has been obtained by 
investigations of the current noise [10], both theoretical [TT ] IT ^ IT3 l lT4] and 
experimental [T51ll6l[T7] . Recently the theoretical interest has turned to the 
full distribution of current fluctuations [18], the Full Counting Statistics 
(FCS) of coherent MAR-transport [T9 ] l20" p T] ■ 

In several types of superconducting junctions, e.g. long diffusive superconducting- 
normal-superconducting (SNS) junctions, quasiparticles lose phase coher- 
ence when traversing the normal part of the junction and both the ac and 
dc Josephson currents are suppressed. This is the regime of incoherent MAR- 
transport. The properties of the current in the incoherent MAR- regime have 
been thoroughly studied, see e.g. |22] and references therein. In addition, 
in several works the current noise has been investigated, both theoretically 
[23 ] l24 ] l2"5] and experimentally [26,27,28,29 . Recently experimental data on 
the third current cumulant in incoherent diffusive SNS junctions was pre- 
sented by Reulet et al [30]. Moreover, Pilgram and the author [31] presented 
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a general theory for FCS of incoherent MAR. The theoretical result for a 
diffusive SNS junction showed qualitative similarities with the experimental 
data [30] . The results in Ref . [31] were mainly derived with a stochastic path 
integral technique |32U33j . however the extension of the FCS for diffusive 
SNS to arbitrary voltage bias was based on a more intuitive picture of inco- 
herent MAR-transport. Appealing to this intuitive picture, in this work we 
first present a detailed derivation and then a thorough analysis of the FCS 
of incoherent MAR in diffusive SNS. Moreover, a formal derivation based 
on the stochastic path integral approach is presented in the appendix. 



2 Model and theory 

We consider a junction with a diffusive normal conductor connected to 
two superconducting terminals. The superconductors have an energy gap A 
and a voltage V is applied between the superconductors. The phase break- 
ing length in the normal conductor is taken to be much shorter than the 
length of the conductor and the transport is consequently incoherent [34l 
135] , We make the simplifying assumptions that the junction is perfectly 
voltage biased and that the scattering is completely elastic. The normal- 
superconducting(NS) interfaces are assumed to be transparent. Due to the 
diffusive nature of the transport, this gives that the effective Andreev re- 
flection is perfect inside the superconducting gap (at quasiparticle energies 
\E\ < A) and negligible outside the gap [3"6]l37]. Throughout the paper we 
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take the temperature to be well below the gap, kT <C A, and consider only 
zero frequency transport properties. 

In the paper we will appeal to a qualitative scattering picture of MAR [3l 
4 , 38J , a more formal derivation based on the approach of Ref . pmi32U35] is 
presented in the appendix. Within the scattering approach, the transport of 
quasiparticles in the junction can be described as MAR-transport in energy 
space |39j : a quasiparticle below the gap injected from the superconductor 
into the normal conductor gain energy by successive Andreev reflections 
before escaping out into the superconductors at energies above the gap. This 
is illustrated to the left in fig. Q] for an applied voltage A < eV < 2 A. Due 




Fig. 1 Left: MAR-transport in energy space: An electron-like quasiparticle is 
injected from the left superconducting terminal at an energy below the gap, E < 
—A. The quasiparticle propagates via MAR to an energy above the gap, E > A, 
at which it escapes out into one of the superconducting terminals. Right: Energy 
windows for the transfer of an effective charge 2e and 3e. 

to the perfect effective Andreev reflection at the normal-superconducting 
interfaces, there are two possible processes by which a quasiparticle, injected 
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from the left, can propagate from filled states below the gap to empty states 
above the gap. 



• for injection energies —A — eV < E < A — 2eV a quasiparticle traverses 
the gap three times (gaining an energy 3eV) as electron — > hole — > electron, 
effectively transferring a charge 3e from left to right. 

• for injection energies A — 2eV < E < — A a quasiparticle traverses the 
gap twice (gaining an energy 2eV) as electron — ► hole, effectively transfer- 
ring a charge 2e from left to right. 

The two transport processes are shown to the right in Fig. [T] Impor- 
tantly, only these processes contribute to the net charge transport. We em- 
phasize that in the scattering processes the net charge transferred across 
the junction is given by the electron charge multiplied by the net number 
of quasiparticle travcrsals or equivalently, the gained energy divided by the 
voltage V. This holds independently on the injection energy and the type 
of quasiparticle injected [IS]. We also note that for quasiparticles injected 
from the right there are, due to the symmetry of the junction, two equivalent 
transport processes. Below, we for simplicity consider only the left injected 
quasiparticles and multiply all the derived transport properties with a factor 
of two. 
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2.1 Generalization to arbitrary voltage 

The scenario presented above can be generalized to arbitrary applied volt- 
age. We have that for a given voltage 

2A 2A 
—-<eV<—, n = 0,1,2,... 1 

n+1 n 

there are two relevant transport processes for left injected quasiparticles. 

• quasiparticles injected at energies 

-A-eV < E < A- (n + l)eV (2) 

traverse the junction n + 2 times, transferring a charge (n + 2)e from left to 
right (gaining an energy (n + 2)eV). 

• quasiparticles injected at energies 

A - (n+ l)eV < E < -A (3) 

traverse the junction n + 1 times, transferring a charge (n + l)e from left 
to right (gaining an energy (n + l)eV). As emphasized above, no other 
processes transferring charge between occupied and unoccupied states are 
possible. Consequently, only these two processes ([2]) and ([3]) contribute to 
the net charge transport. 

2.2 Effective diffusive transport 

The MAR-transport in real space is diffusive. Since the transport is incoher- 
ent, neither the energy of the quasiparticle inside the normal conductor nor 
the electron or hole character of the quasiparticle is of importance for the 
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transport properties. As a consequence, the process where a quasiparticle 
traverses the junction a net n number of times, transferring a net charge 
roe, gives the same (zero frequency) transport statistics as a process where 
a quasiparticle with an effective charge q n — ne is transported through n 
normal diffusive conductors in series. The ro conductors in series have a con- 
ductance G n — G/n, with G the conductance of the individual conductors 
[40] . We point out that this effective model, illustrated in fig. [21 can be seen 
as a straightforward extension and combination of the noise models of refs. 



Fig. 2 Left: MAR-process transporting a charge 3e from left to right. Right: 
Equivalent normal circuit, conductance G3 = G/3, with quasiparticles with effec- 
tive charge 3e. 

The "effective voltage" across the equivalent normal conductor is just 
the energy window available for the process (divided by e). Moreover, each 
of the two transport processes possible for a given voltage are independent. 
Consequently, the the total transport statistics is obtained by just summing 
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up the contributions from the two processes. Importantly, since the FCS, 
or the transport statistics of diffusive conductors is known [42] . the above 
observation allows us to directly derive the FCS of the diffusive supercon- 
ducting junction. 



3 Cumulants of the current 



Considering only positive voltages, we start with the average current. The 
current is the sum of the currents of two transport processes. For a voltage 
2A/(n + 1) < eV < 2A/n we have 



n+2" 



2A 

(n + l)V- — 

e 



G 



n + 1 



(n + l) + 



— -nV 

e 



G 



n + 2 



(n + 2) 



GV 



(4) 



which is just Ohms law. Here, for clarity, we introduced the effective voltages 
for the two processes = V(n + 1) - 2A/e and T4 (2) = 2A/e - nV. As is 
clear from Eq. although the individual processes depend in a nontrivial 
way on the applied voltage via n, the total current is just linear in voltage. 

For the noise, the transport processes are associated with a certain noise 
power P. For a diffusive conductor the noise power for n conductors in series 
is P n = 2eG n /3, where the 1/3 suppression factor comes from the diffusive 
nature of the transport 43 , 44J . In addition, the noise is proportional to the 
effective charge squared. We can thus write the total noise 



s(v) = v^ Pn+1 (*ji) 2 + vw Pn+2 (*^y 
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(n + l)V 



2zA 



G n+ i(n+l) 2 + 



2Z\ 



2A 



nV 



G n+2 {n + 2) 2 



V([n + I] 1 - n[n + 2]) + ([n + 2] - [71 + 1]) 

e 



3 V 2zi/ 



(5) 



which is the known result of Refs. [24,25j. It is interesting to note that again 
the total noise shows no n-dependence, it is a linear function of voltage, with 
a zero voltage offset. However, going to the third cumulant C3, the result 
does depend on n, i.e. it displays a subharmonic gap structure in the form 
of kinks at voltages 2A/em with m = 1,2, ... The expression for the third 
cumulant is given, along the same lines as for current and noise, by 



C 3 (V) = —G [eV{l -n-n 2 ) + A{6 + 4n)] 

1 



In the same way we also get the fourth cumulant 



(6) 



C 4 (V) = — G [eV(-l + 4n + 6n 2 + 2n 3 ) - 2A(7 + 9n + in 2 )] (7) 



The first four cumulants are plotted in fig. [3] We see that, just as for a normal 
diffusive conductor, the first three cumulants are manifestly positive but the 
fourth cumulant is negative. In addition, while the first two cumulants are 
finite at low voltages, the third and fourth cumulants diverge for V — > 0. 
Interestingly, in the experimental data [30] the third cumulant shows a clear 
increase with decreasing voltage, down to a cut-off voltage. In this context 
we note that e.g. inelastic scattering, not accounted for in our model, will 
remove the low voltage divergences, similar to the situation for the noise 

[mum]. 
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Fig. 3 The first four cumulants of the current C P (V) as a function of voltage. 
Note that & = I and C 2 = S. 

The subharmonic gap structure in the voltage dependence of the third 
and fourth cumulant is hardly visible in fig. \3\ However turning to the dif- 
ferential cumulants, the derivatives of the cumulants with respect to the 
voltage, the subharmonic gap structure becomes apparent. The differen- 
tial cumulants are shown in fig. 0) We note that the third and the fourth 
cumulants display a clear step structure. From Eqs. ([6]) and we have 
dC 3 (V) e 3 



dV 15 



G{l-n-n z ) 



(8) 



and 



dd{V) e 4 



■G(-l + 4n + 6n 2 + 2n 3 ) 



(9) 



dV 105 

We point out that in contrast to the normalized cumulants in coherent single 
mode junctions in the tunneling regime [13, 14, 19, 20, 21 , the height of the 
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Fig. 4 The first four differential cumulants dC p (V)/dV as a function of voltage. 
The third and fourth cumulants show a clear step structure. 



steps of the differential cumulants can not directly be related to an effective 
charge. 



4 Full counting statistics 

With the knowledge of the FCS of a normal diffusive conductor [l^], the 
scheme above can readily be extended to all higher cumulants. The cu- 
mulant generating function, from which the cumulants follow by successive 
derivatives with respect to the counting field \i is then given by the sum of 
the generating functions for the two MAR-processes 

F{x) = 2_ G ( (^+l)e^-2^ as . nhVex(n+1) _ 1 
e l \ n + 1 

+ 2A- neV a8inh 2^/ ex{ n+2) _ A ( 10) 

n + 2 / 



12 P. Samuelsson 

where r is the measurement time. The current is / = (e/r)dF / dx\ x= $, the 
noise S = 2(e 2 /T)d 2 F/d X 2 \ x =o, C 3 = (e 3 /T)d 3 F/dx 3 \ x =o etc. We note that 
the distribution of transferred charge Q is given by a Fourier transform 

P{Q) = ! d X e F(lx) - lQx (11) 



Here we however do not further investigate the probability distribution, we 
instead turn to the low and high voltage limits of the cumulant generating 
function. 

5 Voltage limits 

5.1 Low voltage, eV <C A 

In the limit of low voltage, eV <C A, we have n ~ lAjeV and the cumulant 
generating function in Eq. (|10|) is given by 

V 2 G 



F( x ) = T -^-asinhVeX 2Zi / eV - 1 (12) 

This is just the FCS for a normal diffusive conductor with an effective 
charge e(2A/eV) transferred. Importantly, as was found in Ref. [31], this 
result (with different G) holds for an arbitrary incoherent superconducting 
junction. The origin of this universal behavior is that charge transport in 
incoherent superconducting junctions at low voltage bias can be described 
as diffusion in energy space. This energy diffusion was already noted for 
diffusive junctions, in the context of noise, in Refs. [241125], 
We note that the different cumulants C p behave as 

C p oc V 2 ~ p (13) 
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This is different from the coherent (short) diffusive SNS-junctions investi- 
gated in Refs. [45|ri9]. where C p oc V 3/2 - p . 



5.2 High voltage, eV > 2 A 

In the limit of high voltage, eV > 2A, giving n = in Eq. (fT]), the cumulant 
generating function in Eq. (| 10[) is given by 

F( X ) = G ([eV - 2A] asinhVe* -l + A asinhVe 2 * - l) . (14) 

The physical interpretation of this result follows directly from the model: 
the first term is due to single particle transport (x in the exponent) with an 
onset at eV — 2A while the second term is two-particle transport (2x in the 
exponent). The "excess" fluctuations, the difference between the cumulant 
generating functions for the diffusive junction with the terminals in the 
superconducting and normal state respectively, is thus given by 

F exc ( X ) = F(x) - F N (x) = GA (asinhVe 2x - 1 - 2asinhVeX - l) . 

(15) 

We note that this result is independent on voltage. It gives the excess cu- 
mulants I exc = 0, S exc = 4AG/3, C\ xc = 2eAG/5, etc. 

6 Conclusions 

In conclusion, we have investigated the full counting statistics of incoherent 
multiple Andreev reflections in diffusive SNS-junctions. An expression for 
the cumulant generating function has been derived. A careful analysis of the 
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first four cumulants and the low and high voltage limits of the cumulant 
generating function has been presented. The results are valid under the 
assumptions of perfect voltage bias and the absence of inelastic scattering. 
To allow for a quantitative comparison with experiments for the full range 
of applied voltages, a theoretical model where these assumptions can be 
relaxed would be of great interest. 
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8 Appendix 

The cumulant generating function in Eq. f) 1 0[) can be conveniently derived 
in a more formal way using the stochastic path integral approach |32[l33j 
applied to incoherent MAR |31j . We consider for simplicity a wire geometry, 
length L, of the diffusive conductor (for a more general diffusive geometry, 
see Ref. [33 ). The starting point is to discretize the diffusive wire into 
K 1 nodes, contacted via quantum point contacts supporting M 3> 1 
transport modes with transparency T. This discretization is shown in Fig. 
[5j The transparency T and the number of transport modes M are chosen 
to give a conductance of the K nodes in series [MT/K]2e 2 /h equal to the 
conductance G of the wire. 



Title Suppressed Due to Excessive Length 



15 



Q 



j+1,k=K 



Q 



ooocoo oo 



k=1 



k=K 



fj+1,k^-l+1,k 



T j,k Aj,k 



Q 



OQj+1,k=1 



.0 



Pj,k=K 



\ o 

o° e 

o 

oY k= i 



Fig. 5 Left: discretization of the diffusive wire into K nodes. The contact be- 
tween the nodes has transparency T and the contact between the nodes and the 
superconducting electrodes has transparency _T. Right: the node structure along 
a part of the MAR-ladder. The distribution functions fj t k and counting fields Xj.k 
of the nodes labeled with j, k are shown. 



The resulting node structure in energy space, along the MAR-ladder, is 
also shown in Fig. [5] Each node is characterized by a quasiparticlc distribu- 
tion function fj t k and a quasiparticle counting field \j,k, where 1 < k < K 
is the number of the node and j denotes the rung of the MAR-ladder. The 
j:th rung goes between NS-interfaces at energies Ej^i and Ej where Ej = 
E+jeV. We take j = for energies just below the gap — A < E < —A—2eV, 
on the left side. The quasiparticles are thus electrons for rungs with j odd 
and holes for j even. Note also that k is counted along the ladder, upwards 
in energy. 

Following the scheme for FCS of incoherent MAR in Ref. [31] , we write 
the total generating function as a sum (energy integral) over the generating 
functions for different MAR-ladders, which contribute independently to the 
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total generating function. This gives 



/-A oo 
dE 
-A-2eV : ' 



j = -oo 



K-l 



fe=i 



(16) 



where the generating function for a given MAR-ladder is a sum of the gen- 
erating functions S^ k and S^ s along the ladder. Here 

Mr 

S» k = — In [f jtk f jtk+1 + UA1 f jtk+1 ) ^Te i,k+1 ~ i ,k + Ft] 
+ f jtk+1 (l - f jjk ) [Te x **- X w + R] 

+ (1-/^0(1-/^+1)] (17) 

is the generating function |46j for the quantum point contact between nodes 
in the wire, with R = 1 — T. The generating function S^ s describes the 
connection between nodes at the NS-interfaces (471148] . At energies inside 
the gap, \Ej\ < A, we have 

Mt 

S? S = — In [/i.x/j+u + Uk(1 - f j+ i,i) [R A e 2 ^ +x ^- x - K + Rn] 
+ f j+hl (l - f j>K ) [R A e- 2x i +x ^- x ^ +R N ] 

+ (18) 

with Ra = Ra(E, r) the Andreev reflection probability and Rn = Rn(E, r) 
1 — Ra the normal reflection probability, where r is the (normal) trans- 
parency of NS-interface. The counting field Xj is taken equal to \ for j odd 
and zero for j even, i.e. the transferred charge is counted on the right side 
of the junction. It is important to note that \ counts electrical charge while 
Aj : fc counts quasiparticles. 
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Outside the gap, at \Ej\ > A, the Andreev reflection can effectively be 
neglected (as discussed above) and we can for simplicity consider a gener- 
ating function for a perfect normal interface, i.e. a NS-intcrface with the 
superconductor taken in the normal state. This gives 

S™ = ™L In [f. Kfj + f. K{ l - /.) [ e x^, K + i] 
+ W-fj,K)[e-x> +X '- K + 1 ] 

+ _/,.)] (19) 

above the gap, 23^ > Z\ and similarity below the gap. Here fj — f(Ej), 
with f(E) the quasipartice distribution function of the superconducting 
electrodes. For kT <C A considered here, we have f(E) = 1 below the gap 
and f(E) = above the gap. 

To evaluate the generating functions for the ladders we first perform a 
transformation 

Xj,k -> A,-,fc - jx, j = 2, 4, 6... 

A,.a • A,> ;./ j = 1,3,5.... (20) 

[Under the stated conditions, only Aj.fe with j > 1 are of interest]. This 
transformation removes the counting field x from the generating functions 
of the NS-interfaces at energies inside the gap, in Eq. (fT8|) . In addition, the 
generating function for the NS-interface above the gap, Eq. (fTU)) is modified 
as Xj iX- This means that the counting field x IS transferred away from 
the interior of the ladder to the upper boundary. 
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Second, we the note that the structure of the generating functions for the 
quantum point contact between the nodes in the wire, Eq. (|17|). and for the 
NS-interface at energies inside the gap, Eq. (Tl8|) , are formally identical when 
X is transferred away (however with Ra instead of T) . The NS-interface thus 
effectively connects the end- nodes j, K and j + 1, 1 of the two rungs j and 

i + i- 

In addition, we emphasize that the boundary conditions (in the super- 
conductors) for the distribution functions as well as the counting fields [af- 
ter the transformation in Eq. (|20p ] are independent on the electron or hole 
character of the quasiparticles, they only depend on the energy Ej. Con- 
sequently, the fact that distribution functions and counting fields describe 
electrons for j odd and holes for j even becomes irrelevant for the charge 
transfer and one thus only need to consider the position j, k of the node in 
the MAR-ladder. 

Taken together, the MAR ladder can thus be represented as a series of 
nK nodes, with n equal to the number of rungs on the ladder, or equivalently 
the net number of absorbed energy quanta eV when traversing the junction. 
Note that since the resistance of the NS-interface is much smaller than the 
resistance of the diffusive wire, the details of the connection at the NS- 
interface (Ra instead of T) can be neglected. The only difference from a 
discretized normal diffusive wire of length nL is then the boundary condition 
for the counting field at the end of the wire, being nx instead of just x f° r 
a normal wire. 
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Having performed this mapping to a discretized normal wire with renor- 
malized boundary conditions we can continue with the standard procedure 
[33] to derive the total generating function in Eq. (fl6|) . We first denote the 
nodes with a single index k running from 1 to nK (dropping the rung in- 
dex). In the limit if > 1 the difference between distribution functions as 
well as counting fields of two neighboring nodes k and k + 1 is of the order 
1/K and we can expand fk+i = fk + ^fk & n d Afc+i = ^fe + ^^fc- Inserting 
this expansion into the generating function in Eq. (|17p we have to leading 
order in Af% and AXk 

Sj? = ^T[AX k Af k -f k (l-f k )A\l] (21) 

Taking the continuum limit and summing over all generating functions for 
different k we get the generating function for a given energy E as [33J 



s[x > f > x] = ifik I dx 







dx dx ^ ^ \dx) 



(22) 



where we have changed from summation over nodes to an integral over the 
dimensionless position < x < 1 along the wire, i.e / = f(x) and A = A(a;). 

Second, this generating function is varied over all possible configurations 
of / and A. This results in a stochastic path integral formulation for the 
averaged generating function F(x) 

eMF{*x})= fvfV\eMS[ix,f,iX}) (23) 

Under the semiclassical conditions considered here, this path integral can 
be solved in the saddle point approximation. The saddle point equations 
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are given by the functional derivatives 



*S\x,f,X] _ 2 d 
5X dx 



dx 



^=0 
dx 2 



Solving these equations for A and / in terms of Xj inserting the solutions 
back into the generating function and performing the integration over x 
then yields (see Ref. [33] for details) 

t G 

F(x) = asinhVe"* - 1 (25) 
e z n 

describing the FCS of a diffusive wire with conductance G n = G/n and a 
renormalized charge e — > ne. Finally performing the integral over energy in 
Eq. (|16|). taking into account energy window for the two MAR-ladders, gives 
the generating function in Eq. (JTUJ). This concludes the formal derivation. 
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